Anomalous rotational properties of Bose-Einstein condensates in asymmetric traps 
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We study the rotational properties of a Bose-Einstein con- 
densate confined in a rotating liarmonic trap for different trap 
anisotropies. With simple arguments, we obtain the velocity 
field of the quantum fluid for condensates with or without vor- 
tices. While the condensate describes open spiraling trajec- 
tories, on the frame of reference of the rotating trap the fluid 
moves against the trap's rotation. We also find expressions 
for the angular momentum and linear and Thomas-Fermi so- 
lutions for a vortex-less state. In these two limits we find the 
same analytic relation between the shape of the cloud and the 
rotation speed. Our predictions are supported by numerical 
simulations of the mean field Gross-Pitaevskii model. 
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FIG. 1. Flow features of an asymmetric vortex-less con- 
densate with form factor a ~ 2 and rotation speed SI = 1. (a) 
Velocity distribution at t = (b) The rotation of the asym- 
metric condensate defines two circles which define the width 
of the spiraling flow in the following picture, (c) Flow on 
the laboratory frame: the shape of the cloud (dashed line) in 
the Thomas- Fermi limit, and the paths of two test particles 
- inner spirals for 6 rotation and outer spiral for 20 rotation 
periods (1 period is 2Tv/n). (d) Flow on the rotating frame 
using the same case as in (c) but for only one period. The 
motion of the particles is opposite to the trap rotation. 

The question of whether an atomic gas may be a super- 



fluid and the quest for distinct superfluid properties are 
among the central goals of recent research with atomic 
Bose-Einstein condensates (BEG). These targets have 
been largely inspired by previous achievements on "^He 
and its superfluid phase He-II, so that from one known 
property of "classical" superfluids some analog is found 
for the "new" condensates. To name a few phenomena 
which have been found this way we can cite the discov- 
ery of anomalous sound velocities |jl|,|| , the discovery of a 
critical velocity for the beginning of viscous damping , 
the study of vortices and their generation |^,|o| and 
the study of rotational properties such as scissors modes 
|]lT| , ^ or moments of inertia [p^ . 

In this context there has been a great interest in achiev- 
ing Bose-Einstein condensation with rotating traps, a 
task which has been completed recently . The mo- 

tivation of this interest is that the ground state of a con- 
densate in a rotating trap can be forced to host one or 
more vortices depending on the angular speed thus 
providing us with persistent currents which are them- 
selves traces of superfluidity. Although most of the the- 
oretical work regarding the generation and stability of 
vortices has focused on isotropic traps [||J^,p, it is only 
using truly anisotropic traps that one may expect the 
condensate to offer a mechanical response to the rotation 
and produce a vortex. This idea is developed throughout 
the paper. 

There are few relevant results related to non-isotropic 
traps. First jl^, the moment of inertia of an inhomoge- 
neous condensate have been studied assuming that the 
gas is an ideal one, i. e. without interactions. Among 
the most striking results is the formula which relates the 
moment of inertia to the asymmetry of the condensate 
given by the expectation value (x^ — y^). In Ref. ^ the 
dynamics and stability of a vortex in this type of traps 
is studied. And finally in Refs. the ground state of 
a condensate in a rotating asymmetric trap is found nu- 
merically for particular values of the rotation speed and 
their stability is studied. 

In this paper we study the rotational properties of 
a condensate in an rotating anisotropic trap, with and 
without a vortex. We consider the problem within the 
mean field approximation, thus taking interactions into 
account (an essential difference with previous work |Q ) . 
Our treatment is both analytic and numeric: we obtain 
expressions for the most relevant observables, which are 
exact in some limits, and verify them with complex nu- 
merical simulations over a large range of parameters. The 



1 



main predictions of the paper are the explicit expressions 
for the moments of inertia and a characterization of the 
flow of the quantum fluid, whose main features are re- 
flected in Fig. 1. 

The model- We wiU study a two-dimensional conden- 
sate in a harmonic trap whose axes rotate at a certain 
angular speed, fl. The potential may be described by 



(5) 



V (x, t) - {e-'^'"'' ^) =V„{U{-nt)x), 



(1) 



where is the generator of the rotations in the 3?^ space, 
U{9) is an orthogonal transformation which rotates a 
vector an angle around the origin, and Vo(?'i,r2) = 
\{ujlrl+ujlrl). 

We will work on two different frames of reference: the 
laboratory frame {5*, x} which is stationary and the ro- 
tating frame {S,y — C/(— rif)x} which moves with the 
trap. On the flrst frame the zero temperature mean fleld 
theory reads 



idtil) (x, t) 



-]^l^^ + V{x,t)+g\^\^ 



^(x,t). (2) 



Here g characterizes the interaction and is defined in 
terms of the ground state scattering length. It is im- 
portant to remark that our solutions have a well defined 
norm that we may take as equal to the number of parti- 
cles Q, IIV'IP = II0IP = N. To simplify the treatment 
we have assumed a set of units in which h = m = 1. 

We will define a second wave function, this time on S, 
given by (j>{r,t) — tp{U{flt)r,t) with r — (ri,r2). The 
evolution of this function is ruled by 



idtct) - 



-^Ar + Voir)+g\cl)f -ni. 
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Here = —i {rid24> — '''2di<t>) is a representation of the 
angular momentum operator along the z-axis. 

Stationary states. - Let us write Eq. (H) in the modulus- 
phase representation in the rotating frame of reference. 
Defining (0 = v^^*^) obtain the continuity equation 
for the density 



dtp ~ nvp ■ (i?,r) = -V • [H 



(4) 



where v = Im(0V0) = V arg — V8 is almost the 
velocity field of the quantum fiuid in the stationary 
frame. Actually the velocity frame on S is given by 

V = u (~nt)w{u {-vtt)). 

Eq. (pi) has a set of stationary states, which are solu- 
tions of the type (t>{v,t) = e~*'^*(/)(r, 0). These states rep- 
resent configurations of the condensed cloud which main- 
tain their shape and move rigidly with the trap (keep in 
mind that r corresponds to the rotating system). 

In our paper we will be interested in the lowest energy 
state, the ground state. The main assumption through- 
out the paper will be that the lines of constant density 
of a ground state are ellipses of an unknown shape, i. e.. 



which is exact for the ground state in certain limits 
to be discussed later. We will also need a normalized 
anisotropy factor, k = (1 — a)/{\ -\- a) € [—1, -1-1]. 
Using Eqs (|^) and (||), one finds 

h [u) (1 - a) nra + ar^di + rsc^a] O = AG, (6) 

where h = dh\ p / d{v?) . Let us look for solutions corre- 
sponding to an incompressible flow, i.e. AG = 0. Then 

n [(1 - a) rir2] + [aridi + rs^s] = 0, (7) 

Writing the solution of Eq. (^ in the form 

0(ri,r2,<) = -pt + 17 " \ r1r2 + Qvort{ri,r2), (8) 
a + 1 

it is clear that there is still an undetermined part of the 
phase, Qvort which satisfies AO„ort = and VOt,ort -L 
Vp. This phase varies along elliptic paths surrounding 
the origin and carries vorticity in case there is any. 

Once we have found the phase, the density may be 
obtained by solving 



M = 



avp (vof 

2^p 2 



+ Va+gp + i 



(9) 



This can be done analytically both in the non-interacting 
limit, gN — > and on the Thomas- Fermi limit, gN — > 00 
as we will show later. 

States without vorticity.- When the EEC ground state 
has no vorticity we can make Qvort = in Eq. (@) 



O = -pt 



it — —rir2. 

a-\-l 



(10) 



As it was already expected, in the radially symmetric 
case (a = 0) the solution is also radially symmetric, the 
phase is uniform and the velocity field becomes zero. 

For an asymmetric trap the velocity field is different 
from zero but still irrotacional, V x V = 0. The defini- 
tion of V above provides us with analytic expressions for 
the fiow on the laboratory. This flow is a product of two 
rotations around the origin: a circular rotation with an- 
gular speed and an elliptic rotation with a smaller an- 
gular speed ui = — fiVl — fc^ on the opposite sense! This 
is a striking two- fold result. First, on the trap frame S 
the fluid is moving along closed elliptic contours oppo- 
site to the rotation of the trap [See Fig. 1(d)] but with 
a smaller velocity. And second, the composition of both 
movements on the laboratory gives a flow made of slow 
and typically open spiraling paths [See Fig. 1(c)]. 

Roughly speaking, the fact that the actual fiow on the 
laboratory frame is much slower is a signature that the 
inertia of a condensate is smaller than the inertia of a 
classical fiuid. We will prove this calculating the angu- 
lar momentum of the cloud, which is given by a simple 
expression 
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{L.), round ^ j P{r^ Ve), = ^ ^'^1.2^11 ■ (H) 

Here we have used a — (^r\) / Ir"^) ■ As it was expected, 
the moment of inertia, / = {Lz) becomes zero for a 
symmetric trap. Furthermore, it is always smaller than 
the classical value 7 = (r^ + rf). We must also remark 
that that Eq. (jl^) differs from the zero temperature limit 
of Ref. because our method is not perturbative and 
allows the shape of the cloud to depend on the rotation 
speed. 

Inserting Eq. (^0|) into Eq. (^) we find the equation 
of a nonlinear harmonic oscillator with a pair of effective 
frequencies which depend on fi, 

u;l,ff^ujl + n'k{k-2), (12a) 
ujl^ff^ul + Q?k{k + 2), (12b) 

We can recover the asymmetry of the cloud from the ex- 
act solutions of this oscillator on two different limits, the 
non-interacting or linear limit, g ^ 0, and the Thomas- 
Fermi limit. Both expressions are 

OL = a (g -> 0) = uj^,,eff/ujy^eff, (13a) 

aTF = a (5 -> OO) = (W:r,e//M,e//)^ . (13b) 

It is apparent in Eq. ( p^ that rotation emphasizes the 
anisotropy. In Fig. 2 (c-d) we compare the Thomas- Fermi 
and linear approximations with two different numerical 
experiments with an anisotropic trap. The main conclu- 
sion is that Eq. (p^(b)) works extremely well for medium 
to large N . There is a small error which depends on the 
norm of the solution. This error is related to the fact 
that a Thomas-Fermi approximation wipes out any de- 
pendence of the shape of the cloud on the the number of 
particles, N. 

Nucleation of the first vortex.- As the rotation speed is 
increased the fluid adapts by evolving to different states, 
which may involve the nucleation of one or more vortices. 
These vortices arise both as zeros of the density and as 
additive contributions, Qvortex, to the phase of the cloud. 
The phase of the vortices suffers a discontinuous change 
around any closed circuit which encloses them. Hence 
vortices satisfy Feynmann's condition for the quantiza- 
tion of the superfluid flow 

^ Ve„ortex • dl = 27rm, m = 0,±l,±2... (14) 

Following the reasoning developed throughout this paper, 
we should now solve a Poisson equation for Qvortex with 
some restrictions on the direction of VQvortex- Instead 
we will follow a more intuitive path which gives accurate 
results. In the framework of our approximation Ovortex 
will correspond to a flow around the lines of constant den- 
sity of the cloud, as it happens in the symmetric case. In 



elliptic polar coordinates, {ri = ucos9,r2 = y/au sinO} 
the gradient of the phase becomes 

W Q yortex ^ g{u) {- sin cos 9) , (15) 

where g(u) is a decreasing function of the radius of each 
ellipse. This expression may be integrated along elliptic 
contours to find the dependence of the vortex phase on 
the elliptic angle 

Qvortex{d) ~mO + TO— ^ \- sin 26. (16) 

2 (a -I- 1) 

We have found the typical dependence for a symmet- 
ric vortex, m0, plus a correction which comes from the 
asymmetry. Using Eq. ( p^ ) we estimate the new value 
of the angular momentum. It has two contributions, one 
coming from the rotation of the cloud and another one 
coming from the vortex, 

, , {rl~rlf 2mJiri)(ri) , , 

(r^-t-rl) (r2-|-r2) 

In the symmetric trap case the introduction of a vortex 
implies a fixed change of the total angular momentum 

l\Lz = (Lz) cortex - {Lz) ground = ^- The asymmetry 
radically changes this picture and the more asymmetric 
the trap is the smaller the amount angular momentum 
kept in the vortex - a quantity which eventually becomes 
zero in the \k\ — > 1 limit. 




FIG. 2. (a-b) Evolution of the moment of inertia per par- 
ticle — Lz/NQ as a function of the rotation speed and A'^ 
for (a) a radially symmetric trap and (b) an asymmetric trap 
with cj| = 2lji. (c) Moment of inertia versus the angular 
speed of the asymmetric trap (solid line), estimation for a 
vortex-less state (dotted line) and estimation for a state with 
a single vortex (dashed line), (d) Asymmetry of the cloud, a, 
in the asymmetric trap as theoretical prediction (dotted line) , 
numerical result for N — 5 (upper solid line) and A'^ = 175 
(lower solid line). All figures are in adimensional units. 
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Numerical simulations.- Up to this point we have 
obtained several predictions concerning the stationary 
states of a BEC in a rotating trap. Although there are 
no experimental results on this yet, we can contrast our 
theoretical work with numerical simulations of the Gross- 
Pitaevskii equation (GPE). 

First let us recall that the stationary solutions of Eq. 
(0) are critical points of the functional B 



E\d>\ = 



r!L, U, (18) 



subject to the constraint of a fixed norm 
Some of the critical points of Eq. (|l8|) correspond to 
ground states. These states represent the typical con- 
figuration of the gas when it is Bose-condensed in the 
rotating trap. To find the ground states we have mini- 



mized Eq. (18) using the Sobolev's gradients method in 
a discrete Fourier basis. This was performed for a two- 
dimensional radially symmetric trap, uji — uj2, and for an 
asymmetric one, uji — \/2'.^2, varying the only relevant 
parameters, f2 and N , over a wide range. As a result we 
obtained two maps which show the angular momentum 
of the ground state as a function of (fi, iV), and which 
are plotted in Fig. 2(a-b). 

In Fig. 2(a) we recover the discontinuous distribution 
which was found in Rcf . Q . Here the system starts with 
zero angular momentum and remains still until a certain 
angular speed, ili(iV). Beyond ^li{N) a vortex grows 
and we have another plateau on which the angular mo- 
mentum remains constant, = TV. Once more there is 
a critical frequency, ^2{N)^ which marks the nucleation 
of another vortex. From then on the evolution of is a 
piece wise differentiable one. There are still more jumps 
due to the creation of more vortices, but Lz is no longer 
constant in each interval, since the vortices may move to 
accommodate more angular momentum. 

In Fig. 2(b) the picture is quite different. First the 
angular momentum has a non-trivial dependence on the 
number of particles and on the angular speed. It is nei- 
ther zero for a vortex-less state, nor constant for a state 
with one vortex. Indeed, the dependence with respect 
to O is extremely well reproduced by Eq. (|l^) and by 
Eq. (|l7|), which suggests that these are exact laws which 
could be derived by some other mean. 

Conclusions.- We have achieved several goals in this 
work. First, we have shown how to obtain analytic results 
about the macroscopic properties of a rotating conden- 
sate, which leads to unexpected and intuitively appealing 
results for the flow of the condensate. And second, we 
have derived precise laws which relate the anisotropy of 
the condensate, its moment of inertia and the two rele- 
vant parameters of the problem: the norm, A'^, and the 
trap angular speed, Q.. Our predictions have been numer- 
ically confirmed generating maps of L^, (x^) and {y^) as 
a function of (O, A^). 



We hope that these methods will support further work, 
such as the study of linear stability of vortices in asym- 
metric condensates, the development of better approxi- 
mations to the shape of the cloud or even a variational 
analysis which describes the transition between different 
vorticities. Furthermore, our predictions are easily ex- 
tensible to three-dimensional condensates giving many 
results (condensate shape, moment of inertia or energy 
release, for instance) which could be tested in current 
experiments Jl2| . For instance, the discontinuity on the 
asymmetry of the cloud when a vortex is nucleated could 
be used as a non- destructive mean to ensure the existence 
of a phase singularity in the atomic cloud. 

This work has been partially supported by CICYT un- 
der grant PB96-0534. 
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